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We present a general geometrical approach to the problem of escape from a metastable state in 
the presence of noise. The accompanying analysis leads to a simple condition, based on the norm 
of the drift field, for determining whether caustic singularities alter the escape trajectories when 
detailed balance is absent. We apply our methods to systems lacking detailed balance, including a 
nanomagnet with a biaxial magnetic anisotropy and subject to a spin transfer torque. The approach 
described within allows determination of the regions of experimental parameter space that admit 
caustics. 
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Introduction. A noisy dynamical system will occasion¬ 
ally experience large fluctuations that can dramatically 
alter its state. These fluctuations are responsible for a 
wide variety of interesting behaviors, including stochastic 
resonance [1, 2], transport via Brownian ratchets [3, 4], 
logarithmic susceptibility in driven non-adiabatic sys¬ 
tems [5], and Brownian vortices [6, 7]. In the limit of 
weak noise, the system’s dynamical response is deter¬ 
mined by its optimal paths, i.e., the paths along which it 
moves with maximum probability [8, 9]. 

When the system’s stationary probability distribution 
lacks detailed balance, the optimal paths can exhibit un¬ 
usual behavior. In particular, singularities known as 
caustics can develop in the action of the quasistation¬ 
ary density [4, 10-13] resulting from folds and cusps in 
the projection of the Lagrangian manifold of escape tra¬ 
jectories onto the original space of the dynamical vari¬ 
ables [10]. Such singularities cannot occur in the pres¬ 
ence of detailed balance, but in its absence their pres¬ 
ence can significantly alter the behavior of noise-induced 
escape from a static metastable state [10, 12] or a limit 
cycle [14, 15]. Optimal trajectories avoid caustics [10], so 
appearance of a caustic in the vicinity of a stable or sad¬ 
dle point can dramatically alter the escape behavior [12]. 

In this Letter, we reformulate the escape problem and 
in so doing determine conditions under which singulari¬ 
ties in the optimal escape trajectories can appear. This 
leads to a new approach, based on a simple feature of 
the deterministic (i.e., zero-noise) dynamics, toward de¬ 
termining the presence and behavior of caustics. We 
apply our approach, first to a previously studied sys¬ 
tem [12] in which caustics are known to dramatically af¬ 
fect the escape dynamics, and then to a system not stud¬ 
ied from this perspective, namely magnetic reversal in a 
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biaxial nanomagnet subject to both thermal noise and 
spin transfer torque. For the latter system we determine 
the experimental parameter ranges for which caustics are 
likely (and unlikely) to occur. 

Escape in a noisy dynamical system. Consider an 
overdamped particle with position vector x(t) in an n- 
dimensional space. If the particle is subject to both de¬ 
terministic and random forces, its time evolution in the 
general case is described by the Langevin equation 

X = F(x) + v^H • W, (1) 


where F(x) denotes the drift field, W represents a white 
noise process, and the tensor H(x) and scalar e char¬ 
acterize the noise anisotropy and strength, respectively. 
We take e to be sufficiently small so the timescale of a 
successful escape is much longer than that of a single 
excursion from the stable point. 

Diagonalizing (1), we let G be the diagonal matrix 
with Gii = g‘1 corresponding to H; the associated drift 
field vector then has components fi = Fi/gf. The quasi¬ 
stationary distribution p{x.) oc exp[— S'(x)/2e] is given to 
leading order in the zero-noise limit [8] by the solution 
of the variational problem S = min dt, with 

the Friedlin-Wentzell (FW) Lagrangian and associated 
Hamiltonian given by 


Lpw 

Hpw 


/x-F).G-i.(x-F) 

I [(p + f).G.(p + f)-f.G.f 


( 2 ) 

(3) 


Where p is the FW canonical momentum with ele¬ 
ments Pi = dx^Lpw- 

Suppose now that the drift field contains one or more 
stable fixed points, and that the system’s initial state lies 
within the basin of attraction of one of these. It can be 
seen by inspection of (3) that the fluctuational dynamics 
can be mapped onto the problem of a particle with unit 
electrical charge moving on a Riemannian manifold with 
(positive definite) metric tensor G, under the combined 
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FIG. 1: Diagram of momentum ellipse parametrized by 7 . 70 
and TT + 70 correspond to instanton and anti-instanton solu¬ 
tions respectively. 


influence of a magnetic vector potential A = — f [16] 
and electric scalar potential 0 = —f • G • f = “|f|^ [1"^]- 
Hamilton’s equations of motion are then: 

i = G.p + F = G.(p + f) (4) 

P = —• (5) 


The optimal escape paths are zero energy trajectories, 
that is, those along which the Hamiltonian (3) van¬ 
ishes (see, for example, [18]). For any state vector x, 
the set of momenta satisfying the zero-energy condition 
Hpw = 0 can be seen to describe an n-dimensional el¬ 
lipse in momentum space centered at p = —f and with 

axes ai = . 

The momentum ellipse. For simplicity we now confine 
our considerations to two dimensions; the procedure pre¬ 
sented below is straightforwardly generalized to higher 
dimensions by parametrizing an n-dimensional ellipse. 
The 2D momentum ellipse can be parametrized as: 

Pxh) = cos 7-/3; (6) 

9x 

Pvh) = (7) 

9y 

This momentum ellipse defines all possible least-action 
motions accessible to the particle at x traveling along 
an escape trajectory. The usual anti-instanton trajec¬ 
tories (p = 0 ) characterizing zero-noise drift correspond 
to 7 = 7 o = — arctan( instanton trajectories for 

systems satisfying detailed balance correspond to motion 
antiparallel to the drift fields {p = —2/) [ 8 , 18, 19], which 
here corresponds to 7 = tt + 70 . Fig. 1 shows a typical 
diagram of such a momentum ellipse. 

Because the angle 7 parametrizes momentum at each 
point in configuration space, it is convenient to express 
the equations of motion of the particle trajectories solely 
as a function of x and 7. Substituting (6) and (7) into 
the Hamiltonian equations of motion (4) and (5) , we find 


X = fla;|f|GCOS7 

y = fly|f|Gsin7- (8) 

The role of 7 in characterizing the direction of escape is 
apparent. The slope of the escape trajectory is found by 
dividing the second of Eqs. ( 8 ) by the first to yield 

dy/dx = {gylgx) tan7. (9) 

The nature of fluctuational trajeetories. Eqs. ( 8 ) show 
that thermally driven dynamics evolve at the rate of 
the deterministic dynamics rescaled by the noise-induced 
metric, i.e. = |f||, = |F||,_^. The two rates 

become trivially identical for isotropic noise (G = 1) 
where detailed balance is satisfied, because the instanton 
X = —f and anti-instanton x = f trajectories are simply 
sign-reversed. It is somewhat surprising, however, to see 
that it holds more generally. This result can alternatively 
be derived by writing down the effective Lorentz dynam¬ 
ics for a charged particle traveling in both electric and 
magnetic fields 


x = V|f|2-xx (Vxf) (10) 


and noting that upon multiplying by x one obtains 
9t|xp = 9t|fp, again implying that the dynamical speed 
of a particle moving under the influence of noise is equal 
to the norm of the zero-noise drift field. This notion 
has been employed in the literature [ 20 ] to construct 
an efficient numerical scheme (the gMAM method) ca¬ 
pable of computing transition pathways via geometric 
minimization of the FW action, improving on the older 
String method [ 21 , 22 ]. 

Using these results, the Friedlin-Wentzel Lagrangian 
can be rewritten as: 


Lpw 



fP- 

Ig-1 



|F||,_, -2x-G-Vf 


Ifll 

|f||, (1 -cos5') 

( 11 ) 


where we have employed the identities F = G • f and 
X • G“^ • X = |f ||,, and defined 


^ = arccos j ( 1 ^) 

as the angle between the instantaneous escape velocity 
X and the deterministic drift field f at x. Eq. (12) 
shows that, as long as T does not vary too much over 
the course of the escape trajectory, the effective action 
-5(x) = Lpwdt will be dominated by the behavior 
of |f|Q, implying that the norm of the drift field alone 
captures much of the structure of the system’s action. 
This is the central result of the paper. 
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To characterize better the relative importance of |f | vs. 
T on the escape dynamics, consider a closed fluctuational 
path C that moves from a stable fixed point xs to some 
other point xq within its attractive basin, subsequently 
returning (along the anti-instant on trajectory x = f) to 
the fixed point. Using to denote the surface enclosed 
by C, we find 


S'c = ^ (|f 1 ^ - X ■ f) dt = ^ ds|f Ig - J df 2 c • (V x f) . 

(13) 

In the second equality we changed integration variables 
from t to s, where s denotes distance traveled along the 
trajectory, and used the earlier result that |x| = |f| for 
all systems. We also used Stokes’ theorem to rewrite the 
effect of the geometrical phase due to T in terms of its 
equivalent surface integral. This equation shows how the 
presence of a nongradient field (which possesses nonzero 
curl) contributes to the action, and therefore helps deter¬ 
mine the trajectory, of an optimal path. While enclosing 
a larger ‘magnetic flux’ f dllc * (V x f) acts to reduce the 
path’s action (where the flux term is positive), doing so 
can also move the particle into regions of configuration 
space with larger |f|, which acts to increase the action. 
The least action trajectory therefore optimizes the rela¬ 
tive balance between these two terms. 

There are two immediate consequences of (13). The 
first is that in gradient systems, where f = —VyJJ for 
some smooth potential function U(x) (so V x f = 0), 
and where the instanton trajectory satisfies x = —f, we 
recover the well-known action (see, for example, [18]) for 


the entire trajectory: Sc 
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U(xo) - U(xs) 


. The sec¬ 


ond, more general consequence is that any sudden change 
in |f|, as system parameters are varied, can dramatically 
alter the balance between the drift norm and the mag¬ 
netic flux, and hence the structure of the optimal paths. 

Applications. It has been observed [ 10 , 12 ] that noisy 
systems with nongradient deterministic dynamics con¬ 
taining tunable parameters can dramatically change their 
escape behaviors, in a way reminiscent of a broken sym¬ 
metry transition, when certain critical parameter thresh¬ 
olds are crossed. Consider, for example, the following 
well-studied system [ 12 ] with a tunable parameter a: 


X = x(l — x^ — ay^) -h y/eWx (14) 

y = -y{l + x^) + CeWy, (15) 


where the additive white noise is isotropic and the drift 
field is nongradient for all a 7 ^ 1. For any a, the system 
has two stable fixed points at (± 1 , 0 ) and one hyperbolic 
(saddle) fixed point at (0, 0). We consider the optimal es¬ 
cape trajectory from ( 1 , 0 ) to ( 0 , 0 ), which flows along the 
x-axis for a < A. Above a = 4, however, caustics appear 
in the basin of attraction of the stable point [ 12 ], focus¬ 
ing to a point on the x-axis between 0 and 1 , as shown 
in Fig. 2 . As a consequence, the optimal escape trajec¬ 
tory bifurcates into two off-axis trajectories. Analysis of 
the norm of the drift field reveals that in transitioning 
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FIG. 2: Contour plot of the norm of the drift field from [12], 
for 0 = 3 (left) and o = 5 (right). Optimal escape trajectories 
are shown in red, and other instanton trajectories in green. 
For 0 = 3 , instanton trajectories do not cross, and the escape 
trajectory lies along the a:-axis. Correspondingly, two global 
minima in the norm of the drift field are present at the unsta¬ 
ble (0,0) and stable (1,0) equilibria, respectively, along with 
a saddle close to the midpoint of the x-axis. For o = 5 , in¬ 
stanton trajectories cross, indicating the presence of a caustic 
(not shown). Correspondingly, the saddle in the norm of the 
drift field for a = 3 is now a local maximum, with two new 
local minima appearing off the x-axis. Consequently, there 
are now two symmetrical off-axis optimal escape trajectories 
that follow the new off-axis minima in the norm of the drift 
field (and in so doing avoid the caustic). 


across the critical threshold ac = 4, the structure of the 
extrema of f changes abruptly. For a < 4, the norm ex¬ 
hibits two global minima at ( 1 , 0 ) and ( 0 , 0 ), along with a 
saddle near the midpoint of the x-axis. For a > 4, how¬ 
ever, two local minima appear symmetrically displaced 
off the x-axis, with the previous on-axis saddle now a 
local maximum (Fig. 2 ). According to (11), these new 
minima lower the effective action of off-axis escape tra¬ 
jectories, in accordance with observations. Deviation of 
the escape trajectories from the exact minima are due to 
the ‘magnetic flux’ term in (13). 

A more interesting application is to a physically rel¬ 
evant model lacking detailed balance: the stochastic 
Landau-Lifshitz-Gilbert-Slonczewski (sLLGS) equation 
governing the evolution of a unit magnetization vec¬ 
tor subject to a combination of both gradient and non¬ 
gradient torques. This equation reads: 


rfii = Ai{m) -I- o Hth.k, (16) 

where the drift vector A(m) and diffusion matrix B(m) 
are given by 

A(m) = m X heff — am x (m x heff) 

— aim X (m x n^), (17) 

Bik{m) = VC[-€ijkrnj - a{mimk - Sik)] (18) 

and the equation is interpreted in the Stratonovich sense. 
The first term in A(m) corresponds to magnetization 
precession about a local magnetic field, with heff = 
— Vme(m) a conservative vector field. Here e(m) is the 
energy landscape of the magnetic system under study. 
The second term in A(m) is a phenomenological damping 
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term, with the damping constant a typically ^ 0(10“^). 
The third is a nongradient term corresponding to spin- 
angular momentum per unit time injected via a current I 
into the macrospin along an arbitrary polarization direc¬ 
tion hp [23]. Although the diffusion matrix B(m) (with 
C the diffusion constant) appears state-dependent, it can 
be shown (e.g., by rewriting the dynamics in spherical co¬ 
ordinates) to correspond to isotropic, state-independent 
noise. 

In the absence of applied currents (/ = 0), the ffuctu- 
ational trajectories are determined by the energy land¬ 
scape e(m) and do not cross. In the presence of a non¬ 
vanishing current /, however, detailed balance is absent, 
and it therefore becomes important to determine how 
this new feature may — or may not — alter the escape 
dynamics. The methods developed above allow us to an¬ 
alyze this problem by examining the norm of the total 
drift field governing the macrospin dynamics. To lowest 
order in a, it is 

|A(m)p = |mx heffp + 2 (a/rip- (mx heff) + 0(a^). (19) 

We wish to determine under which conditions the con¬ 
servative precessional contribution dominates the non¬ 
gradient contribution. When this occurs, the extrema of 
|A(m)p corresponding to / = 0 do not change signifi¬ 
cantly when / > 0. By the approach developed above, 
this implies that the most probable escape paths should 
not differ significantly from the reverse-drift instanton 
paths of the purely gradient case. 

We consider for simplicity the case of a biaxial 
macrospin subject to the energy landscape 

e(m) = Dml - ml , (20) 

defined on the surface of the unit sphere (|m| = 1 ). 
Here D > 0 is the ratio of hard axis to easy axis 
anisotropies [24, 25]. 

The calculation analyzing the relative contributions of 
the gradient and nongradient terms in (19) is done in 
detail in the Supplementary Notes [26]. We find that 
the presence of caustics is determined by the fixed tilt 
angle uj that rip makes with the z-axis. For any D, the 
escape dynamics are unaffected by caustics for a wide 
energy range in the small tilt case (tancj ^ but 

in principle the large tilt case (tancj ^ a”^) can show 
substantially different behavior. 

We tested these predictions by numerically integrat¬ 
ing the FW dynamics associated with the stochastic pro¬ 
cess (16) for several different cases. Results for D = 0 are 
shown in the Supplementary Notes [26]. The absence of 
caustics, for both zero and small tilt, is apparent. Here 


we present the more interesting case of nonvanishing D. 
Fig. 3 shows results for D = 20, / = 0.8 /c, and three tilts 
cj = 0, 0.1 and 0.2b Oc^ where Oc = arctan(l/V^) [26] 
and Ic is the critical applied current above which the 
entire > 0 region becomes unstable, eliminating the 
bistability of the magnetic system. These tilt values (and 
even larger ones) can be realized in current experiments 
on orthogonal spin-valve devices [27]. The optimal tra¬ 
jectories change as tilt is increased, even though the exit 
point remains essentially unchanged. However, numeri¬ 
cal results for larger tilts, up to the critical tilt, do not 
exhibit crossing of escape trajectories, suggesting that 
caustics have not formed near escape paths. The only 
major difference is the number of precessions the system 
undergoes before reaching the separatrix (i.e., the bound¬ 
ary of the domain of attraction of the stable fixed point, 
which here corresponds to e = 0 [26]). 

To summarize, an analysis of the norm of the drift 
field indicates that for most tilts, caustics do not appear 
within the escape region, and so — perhaps surprisingly 
— the loss of detailed balance does not qualitatively al¬ 
ter the escape dynamics. However, on closer examination 
one finds that for any tilt, there remain regions — where 
rriz ~ 1 or, separately, where e ^ 0 — where this conclu¬ 
sion breaks down, because the precessional contribution 
to the drift norm vanishes at = 1 and at e = 0. That 
is, two regions will always exist in the magnet’s configura¬ 
tion space where the (nongradient) spin transfer torque 
term dominates the (gradient) precessional term. We 
compute the width of these regions in the Supplemen¬ 
tary Notes and show that they are sufficiently small that 
their effect on the escape dynamics is negligible. 

In the case when the tilt is large (tancj ^ (a“^), how¬ 
ever, the nongradient term dominates the gradient term 
in all of configuration space, and so caustics may appear. 
In order to determine whether they do requires a more 
lengthy analysis, in which the fixed point structure of the 
drift field norm must be analyzed to determine whether 
it changes. We defer such an analysis to future work. 
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SI: RELATIVE CONTRIBUTIONS OF GRADIENT AND NONGRADIENT 


TERMS TO THE DRIFT FIELD NORM FOR A BIAXIAL MACROSPIN 

In this section we calculate the relative contributions of the gradient (precessional) and 
nongradient (spin torque) terms in Eq. (19) in the text. In the absence of applied currents 
and damping, the motion is purely precessional. Focusing solely on this precessional motion, 
we can divide the unit magnetization sphere into four disjoint dynamical regions, as shown 
in Fig. 1. When magnetization energy is conserved, the surface of the sphere can be divided 


z 



FIG. 1 : Constant energy contours of the unit magnetization sphere. Red and blue contours corre¬ 
spond to e < 0 and e > 0 precessions respectively. Dashed black lines are the e = 0 separatrices. 


into trajectories with e > 0 and those with e < 0. The e < 0 surface can be further 
subdivided into orbits with > 0 and < 0 , and the e > 0 surface similarly subdivided 
into orbits with nrix > 0 and rrix < 0. These regions are separated by a common separatrix 
rrix = corresponding to a degenerate e = 0 dynamical state whose trajectory lies 

on a plane tilted at an angle 6c = arctan(l/v^) with respect to the y-z plane. 

However, the above description must be modihed when damping and current are taken 
into account, since either will cause the magnetization energy e(m) to change over time. 
Using (16)-(18) from the text and the definition heg = — Vine(m), we find 

e = Vme • in = -a [m X (heff -|- Ihp)] ■ (m x heff) . (1) 

Several qualitative conclusions can be drawn from (1). When / = 0, dissipation relaxes the 
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system to one of the two negative energy basins in an effort to align the magnetization with 
the ;2-axis. The addition of current (which we will conventionally choose to be negative, 
/ < 0) leads to a progressive destabilization of the > 0 negative energy basin in favor 
of that at < 0. More generally, the presence of current alters the boundaries of the 
attractive stable basins, even though a proper energy function that includes the spin-torque 
terms cannot be written. However, one can show that there exists a critical value Dq ~ 5.09 
that separates two distinctly different qualitative regimes [1]. In the D < Dq case, the current 
alters the e = 0 separatrix boundary, while in the D > Dq case, the e = 0 boundary remains 
intact while the stable fixed points are modified by the current’s action. Furthermore, above 
a critical current Iq (to be discussed further below) the entire > 0 region becomes 
unstable, eliminating the bistability of the magnetic system. 

We want to determine under which conditions the nongradient current term can alter the 
fixed point structure of the norm of the drift field. Because of the smallness of a, energy 
changes slowly, so one can still gain useful information by studying the norm of the drift 
along a constant negative energy trajectory. Along each energy level set e one can compute 
the maxima and minima of the two terms appearing on the RHS of (19). This allows us to 
determine the range of experimental parameters for which the precessional term dominates, 
and therefore when effects due to nongradient terms can safely be ignored. Using (20) 
and the definition of he®, we find the minimum of the precessional term along a particular 
constant energy contour to be: 

min{|m x heff|^} = 4min{e(T> — e) -|- (T> -|- l)m^} = 4|e|(l — |e|), (2) 

where we used the fact that min{m^} = —e along a constant negative energy contour. 

We next consider the maximum of the nongradient current term in (19) and, for def¬ 
initeness, we limit ourselves to cases where hp lies in the x-z plane. Let u denote the 
fixed tilt angle that fip makes with the .s-axis. Then for small a the critical currents 
are Iq = (2/7r)-\/D(D'9M)/cos a; when D > Dq and Iq = (i9-|-2)/2cosa; when D < Dq [1]. 
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Using these results, we find 


max{2«/np • (m x heg)} = (da/c cos u))msix{my {Dnix + tana;)} 

4«-\/|e|(l — |e|)tana;, if D = 0. 

^f^D{D + l){D+\e\)\e\Q{u), if D > Dq. 
2a{D + 2)^{D + \e\)\e\Q{u), if D < Do. 


= < 


(3) 


where 


Q{uj) = maxn 


' 1 - 


{D + l)m 


D+\e\ 



-1 + 


tana; 


(4) 


and we set the current equal to Ic in each case to obtain an upper bound. 

The above equations reveal a quantitative difference between the small tilt (tana; <C a~^) 
and large tilt (tana; ^ a~^) cases. Consider first the case D = 0. By (2) and (3), the 
precessional term dominates the nongradient term when e(l — e) ^ tan^ u. For the small 
tilt case, as long as the magnetization energy is not too close to the separatrix (e = 0) or 
the rUz = 1 pole (|e| = 1), the precessional term dominates and no caustics will be present. 
More precisely, the precessional term dominates within the energy range a^tan^o; < |e| < 
1 — tt^tan^o;. When the tilt is larger, on the other hand, the precessional term dominates 
for most parameter values, and the presence of caustics may affect the reversal dynamics. 
Whether they do requires a detailed analysis of the particular case of interest. 

This conclusion remains qualitatively unchanged when D > Dq, with the main quanti¬ 
tative change being that the energy range in which caustics cannot appear in the small tilt 
case is now Da tan a; < |e| < 1 — Da tan a;. For D < Dq, we need to further subdivide into 
the D > |e| and D < |e| cases. The latter smoothly goes to the D = 0 case as D decreases, 
while the former gives the same result as the D > Dq case. Our overall conclusion is that 
for any D the escape dynamics are largely unaffected for a wide energy range in the small 
tilt case, but in principle the large tilt case can show substantially different behavior. 


S2: ANALYTICAL SOLUTION OF THE w = 0 CASE 

In this section we analytically solve for the Hamiltonian dynamics arising from the FW ac¬ 
tion for the case of a; = 0, with arbitrary applied current. In this case, the dynamics be¬ 
come effectively one-dimensional and the dynamical equations can be exactly solved. The 


4 



FW Lagrangian for the thermally activated dynamics of the macrospin, expressed in spher¬ 
ical coordinates (6*,^), is: 


Lpw — - 


2 L 

Passing to the Hamiltonian formalism we have: 

1 




Hpw — 


Pe + 


pI 


sim(0) 


+ P0f0+P<l>U, 


where p 0 and are the conjugate momenta of 9 and 0. 
Hamilton’s equations are then 


0 = 


Pe = 
P4, = 


+ f4> 


Pe + fe 
P4> 

sin^d 

“ Pedefe - P^def^ 

Sin 0 ^ 

-Ped^fe - P4,d4,f4„ 


(5) 

( 6 ) 


(7) 

( 8 ) 

(9) 

( 10 ) 


where so far all results are completely general and extensible to models of any complexity 
(anisotropy ratio D 0, positive tilts between axes, etc.). For a macrospin with both D = 0 
and (u = 0, the drift vector field of the dynamics (expressed in terms of 6*,0) can be obtained 
by rewriting Eqs. (16)-(18) from the text in spherical coordinates. One then obtains 


fo = —a{I + cos 9) sin 0 
= -cos 9. 


(11) 

( 12 ) 


Because (p is cyclic, p^, = c with c constant. For a zero-energy trajectory of (6), po is given 
by: 

^ ' TiF ^F^l (13) 


Pe = -fe 


l±il-c 2^ 


ac 


' fe ' ffesm^9^ 

This solution is physically valid as long as the radicand is greater than or equal to zero for 
all values of 9 in the range of interest. This implies that c = 0. The conjugate coordinate p 
of p^ then evolves along the deterministic drift field during the escape process, i.e., <j) = f^, 
and therefore p 0 = —f 0 [l ± 1]. The = 0 solution corresponds as usual to relaxation to 
equilibrium, and the p 0 = —2/^ solution corresponds to instanton escape, antiparallel to the 
drift, along the 9 coordinate. However, during this process the precessional dynamics of p 
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remain unaltered, so escape from the stable well does not merely correspond to time-reversal 
of the deterministic dynamics. 

Having solved for the general momenta in terms of the spherical coordinates, we can 
analytically derive the relation between (p and 0 along the escape trajectories. Dividing (p 
by 6 from Hamilton’s equations, we hnd: 


m = 


M- 


COS 9 


a 


1 


a{l-P) 


(/ -|- cos 6 ) sin 9 

log [tan( 0 / 2 )] -h / log 


2 {I + cos 9) 
sin 6 * 


(14) 


where the integrand diverges at the 0 = 0 pole due to degeneracy in the cp coordinate, and 
at the 9 = arccos(—/) separatrix where the escape trajectory is nondifferentiable. 

We tested these results and predictions by numerically integrating the FW dynamics 
associated with with the stochastic process (16) for several different cases. The case D 7 ^ 0 
is discussed in the text; here we show results for D = 0 , with both zero and small tilt. 

Fig. 2 shows the optimal escape trajectory for D = 0, a; = 0, and / = 0.3/c, where, as 
noted in the text, Ic is the critical applied current above which the entire > 0 region 
becomes unstable, eliminating the bistability of the magnetic system. As shown above,, the 
FW dynamics are exactly solvable in this case. 



FIG. 2: Escape trajectory for a macrospin model with a — 0.01, D — 0, I — 0.3 Ic and cj = 0. 
Blue line shows the least-action result of numerical integration of the FW dynamics. Dashed red 
line is the analytical result (14). 


We next consider the case of small tilt. In accord with our analysis, the numerical results 
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demonstrate that the escape dynamics are essentially unaltered (Fig. 3). In the limit of 
sufficiently small a, critical currents scale trivially with the tilt (/^ = coscj) [1]. By 
including this correction to the analytical expression for the uj = 0 escape model, we can fit 
the data well, even though the system no longer obeys detailed balance. 



FIG. 3: Escape trajectory for a macrospin model with a = 0.01, D = 0^ I = 0.3 Ic and cj = 0.3 9c 
(defined in text). The purple line shows the least-action result of numerical integration of the 
FW dynamics. The dashed red line is the analytical result (14). 


[1] D. Pinna, A.D. Kent and D.L. Stein, Phys. Rev. B 88, 104405 (2013). 
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